Massive Gravity Coupled to Galileons is Ghost Free 
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It is possible to couple Dirac-Born-Infeld (DBI) scalars possessing generalized Galilean internal 
shift symmetries (Galileons) to nonlinear massive gravity in four dimensions, in such a manner that 
the interactions maintain the Galilean symmetry. Such a construction is of interest because it is not 
possible to couple such fields to massless General Relativity in the same way. We show that this 
theory has the primary constraint necessary to eliminate the Boulware-Deser ghost, thus preserving 
the attractive properties of both the Galileons and ghost-free massive gravity. 



Introduction and Outline: It has long been known 
that the Fierz-Pauli action [1] provides a consistent de- 
scription of the linear fluctuations of a massive graviton 
in flat spacetime. Nonlinear theories of massive gravity 
tend to suffer from an instability known as the Boulware- 
Deser (BD) ghost [2]. A general Lagrangian for nonlin- 
ear massive gravity can be formulated by introducing a 
non-dynamical reference metric g^ (e.g. the Minkowski 
one, g^v = r\^v) and constructing a potential of the form 
V{g^ a g aiy ). The potential explicitly breaks diffeomor- 
phism invariance and it is expected that the theory gener- 
ally propagates 12 phase-space degrees of freedom, rather 
than the 10 necessary to describe a massive graviton. The 
extra degree of freedom is the BD ghost. 

The problem of finding a ghost-free nonlinear theory 
was only recently solved by de Rham, Gabadaze and Tol- 
ley (dRGT) [3, 4]. dRGT theory is a 3-parametcr family 
of potentials whose special structure ensures that there is 
a dynamical constraint which removes the ghost degree 
of freedom. This has been demonstrated by explicitly 
counting degrees of freedom in the Hamiltonian formal- 
ism [5, 6], and through other methods [7-10]. 

The diffcomorphism invariance broken by the mass 
term can be restored through the Stuckelberg method 
[11], which involves introducing four auxiliary scalars 
(j) A {x) through the replacement g^ v — s- d^^d^^r/AB- 
The Stuckelberg fields are pure gauge and the original 
theory is recovered by choosing unitary gauge <f) A = 
S^x^. In this formulation, the scalars can be regarded 
as the embedding mapping of a sigma model £ — > Ai, 
where both £ and M are four dimensional Minkowski 
space. There is an internal Poincare symmetry corre- 
sponding to the isometries of the target space. From this 
point of view, the dynamical metric g^^x) is a world vol- 
ume metric living on S. 
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The target space may be higher dimensional, and need 
not be flat - we may generalize the sigma model to map 
to an arbitrary target space of dimension D > 4 with 
coordinates 4> A (so that now A,B, - ■■ run over D values) 
and a fixed target space metric Gab (</>), 



9^ = d^d^GAB^). 



(1) 



There are now (dim M) — 4 embedding fields which can- 
not be gauged away and these become physical Dirac- 
Born-Infeld (DBI) scalars coupled to the physical metric 
through the dRGT potential [12]. Apart from the dRGT 
terms, curvature invariants constructed solely from g^ 
and extrinsic curvatures of the embedding can be in- 
cluded in the action. The leading term in the derivative 
expansion is the DBI action ~ J c\ A XyJ—g, and higher 
Lovelock invariants give Galileons [13, 14]. The theory 
will possess a Galileon-like internal symmetry for every 
isometry of Gab, an d the resulting Galileons will be the 
generalized curved space Galileons discussed in [15-18]. 

Apart from generalizing dRGT, the construction of [12] 
is of interest because it provides a method of coupling 
the Galileons to (massive) gravity while preserving the 
Galileon invariance. When coupling to ordinary massless 
gravity, non-minimal couplings must be added to ensure 
second-order equations of motion, and the Galileon sym- 
metry is broken [19, 20]. In the present construction, 
there is no such problem, suggesting that the Galileons 
more naturally couple to a massive graviton. 

In [12], it was shown that the theory is ghost- free, for a 
flat target space metric, in the decoupling limit, and for 
a certain simplifying choice of parameters. In this letter, 
using methods similar to those of [10], we demonstrate 
that the full theory, for any target space metric Gab, to 
all orders beyond the decoupling limit, and for all choices 
of parameters, has the primary constraint necessary to 
eliminate the Boulware-Deser ghost. 

The model: The dynamical variables are the physical 
metric g^ and the D scalars (jy^, which appear through 
the induced metric (1). The action is 



S' = <Seh[3] +Smix[g,g] +S'Gameon[£?] 



(2) 
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Here 5eh[5] is the Einstein- Hilbert action for g^ Vl with 
a possible cosmological constant A, 

hf 2 r 

Seh [9] = ^fj^x^—g (R[g] - 2A) . (3) 
The action mixing the two metrics is 

S mi *[9, g] = — \- E (v/F 1 ?) , (4) 

n— 1 

where y <7 _1 5 is the matrix square root of the matrix 
Q^Qavi an d S n {M) of a matrix M are the symmetric 
polynomials 1 S n (M) = M [ % ■ ■ ■ M''" p ] „ . The /?„ are 
three free parameters (one combination of which is re- 
dundant with the mass to). iSGaiiieon [g] stands for any 
Lagrangian constructed from diffcomorphism invariants 
of g (and extrinsic curvatures of the embedding) whose 
equations of motion remain second order in time deriva- 
tives. The possible terms in if>Gaiiieon [g] are the Love- 
lock invariants and their boundary terms (see [13] and 
Sec. IV. B of [14] for a discussion). The structure of the 
dRGT-DBI coupled system (2) is nearly identical to that 
of ghost-free bi-gravity [21], the difference being that one 
of the two metrics is induced from a target space, and so 
it fundamentally depends on the embedding scalars. 

Following [10], we will find it convenient to write the 
theory in vierbein form 2 . We write the physical metric 
and induced metric in terms of vierbeins E A = E^dx 11 , 
E A = E A dx», 

g ia ,=E li A E v B nAB, fa = E A E„ B r, AB , (5) 

where t\ab is the 4 dimensional Minkowski metric. For 
the induced metric g^, we write the vierbein in an upper 
triangular form 

v - (: . w 

Here N and N 1 are ADM lapse and shift variables, and 
ej a is an upper triangular spatial dreibein for the spatial 
part of the induced metric and e 1 a its inverse transpose 
(in what follows are spatial coordinate indices 

raised and lowered with the spatial metric g~ij, and a, b, . . . 
are spatial Lorentz indices raised and lowered with 5 a b). 
These are obtained in terms of <j) A by solving 

goo = WGabW = -N 2 + N*Ni 
g i = j> A di4> B G AB (<t>)=Ni 
ga = di (j) A dj(f)® Gab (4 1 ) - e, a e 3 b 5 ab . (7) 



1 Our anti-symmetrization weight is . . . fj, n ] = • • 'Mn + 
• ••). Sec appendix A of [10] for more details on the symmetric 
polynomials. 

2 See also [22, 23] for covariant methods of degree of freedom count- 
ing in the vierbein formulation of massive gravity. 



The upper triangular vierbein (6) has 10 components, 
and is merely a repackaging of the 10 components of 
g~p.v (which in turn depend on the <fi A ). 

For the physical metric g^, we parameterize its 16 
component vierbein as a local Lorentz transformation 
(LLT) A, which has 6 components, times a vierbein E 
which is constrained in some way so that it has only 10 
components, 

£/ = K A B E, B (8) 

The freedom to choose the constraints for E allows us 
to make different aspects of the theory manifest. The 
mixing term (4), in terms of vierbeins, takes the form 

M 2 p m 2 A p n (w) 
^ mix " 8 ^n!(4-n)! mix ' 

S £l = Je ABCD E A AE B AE c AE D , 

4t = / e A BCDE A AE B AE c AE D , 

= / cabcdE a AE b AE c AE d . (9) 

The dynamical vierbein has 16 components, 6 more 
than the metric. If wc choose the 6 constraints which E 
must satisfy to be the symmetry condition, 

E^aE^ = 0, (10) 

then we can show using the arguments in [10] (see also 
[24] for subtleties) that the extra 6 fields in A are aux- 
iliary fields which are eliminated by their own equations 
of motion, setting A = 1, and the resulting theory is 
dynamically equivalent to the metric formulation (2). 
Instead, we take E to be of upper triangular form 

V - (: ^) ■ <n> 

Here the spatial dreibein is arbitrary, containing 9 
components. The LLT A in (8) depends now on 3 boost 
parameters p a and can be written as 

where 7 = ^/l + p a p a . Using this decomposition, the 16 
component vierbein E^ A is parameterized in terms of the 
3 components of p a , one N, 3 components of N l and the 
9 components of e, '. 

Hamiltonian analysis: We start the Hamiltonian 
analysis by Legendre transforming with respect to the 
spatial vierbein e i °, defining canonical momenta 7r l a = 
jjgTr. Since S'mix contains no time derivatives of the 
physical metric, and ^Gaiiioon has no dependence on the 
physical metric at all, the expressions for the canonical 
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momenta are the same as their GR counterparts. In par- 
ticular, there will be three primary constraints 



Tab = e.i[ a ^\] = • 



(13) 



In GR, these are first class constraints which generate 
local rotations. 

The Einstein- Hilbert part of the action takes the form 3 

Seh = / d 4 ^ 7r* e 4 ° - \ X ab V ab - NC (e, tt) - N% (e, tt) 

(14) 

The anti-symmetric A ab holds the three Lagrange mul- 
tipliers for the three primary constraints (13). The N 
and N l appear as Lagrange multipliers enforcing respec- 
tively the Hamiltonian and momentum constraints of 
GR: C = 0, C t = 0. 

We now look at the mixing terms (9). 
The contributions to £ m ix arc of the form 



<E P C E C 



D 



containing various 



numbers of copies of E and E. From (6), (11) and (12), 
we see that the [i = components of E^ and E^ A 
are strictly linear in their respective lapses and shifts 
and the fj, = i components are independent of the lapse 
and shift. Therefore, due to the anti-symmetry of the 
cpsilons, the entire mixing term is linear in the lapses 
and shifts, so we may write 

£mix = -iVC mix (e, e,p) - iV l C m i X ,i(e, e,p) - 7VC mix (e, e,p) 
- iV 4 C miXji (e, e,p) - U m ^(e, e,p) . (15) 

The lapse and shift remain as Lagrange multipliers, 
and the p a appear algebraically. We now solve the con- 
straint enforced by N l for the p a : Ci+C m i X) i = => p a = 
p a (e, e, tt). Plugging back into the action we obtain 



S = I d 4 xTr l a e. l a 



1 



A Qb P Qb -iV[C(e,^+C mix (e,e,7r)] 



- iVC mix (e, e, tt) - iV 4 C miXii (e, e, tt) 

^mi X (^, 6, Tt) -\- ^Galileon 

(e,N,N*) 



(16) 



It remains to Legcndrc transform with respect to the 
scalars (f> A , which appear through the dependence of N, 
N l and e^, as determined by (7). In order to avoid deal- 
ing with the complications of diffeomorphism invariancc, 
we first fix unitary gauge, setting the first four fields equal 
to the space-time coordinates: M = x M (this can be done 
consistently in the action, since the missing equations of 
motion are implied by the remaining equations). The La- 
grangian (16) then depends on the remaining (dim Ai)— 4 
scalars and their time derivatives. Crucially, we see from 



3 See [25] or Appendix B of [10] for details of the Hamiltonian 
formulation of GR in vierbein form. 



(7) that while N and iVj depend on time derivatives of 
the scalars, the (ii a, s do not, and this in turn implies that 
the momenta conjugate to the scalars are independent of 
the dynamical lapse N. Thus, when the scalar velocities 
are eliminated in terms of the momenta, the action will 
remain linear in N. (If this were not the case, the lapse 
would no longer be a Lagrange multiplier, but would in- 
stead become an auxiliary field which does not impose a 
constraint on the remaining variables.) The phase space 
is spanned by the 9 independent components of e°~, the 
physical scalars, and the canonical momenta. Since the 
interaction terms break the local rotation invariance of 
GR, the 3 primary constraints (13) associated with the 
local rotations will generate secondary constraints and 
form 3 second class pairs, thus removing 3 degrees of free- 
dom. The constraint enforced by N is precisely the spe- 
cial primary constraint needed to remove the Boulware- 
Deser sixth degree of freedom, leaving 5 degrees of free- 
dom for the massive graviton. Analogously to what hap- 
pens in massive gravity, we expect this special primary 
constraint to generate a secondary constraint to eliminate 
the ghost's conjugate momentum [6]. 

We have implicitly assumed that ^Gaiiieon can be writ- 
ten in such a way that the (dim.M) — 4 unitary gauge 
scalar fields appear with at most first time derivatives, so 
that we may define canonical momenta in the usual way. 
This is not immediately obvious, because the higher order 
Galileons in S'Gaiiieon possess higher derivative interac- 
tions. However, the higher derivative interactions within 
^Gaiiieon are special in that they generate equations of 
motion which are no higher than second order in time. 
This means it should be possible, after integrations by 
parts, to express these Lagrangians in terms of first time 
derivates only (though we shouldn't expect to be able to 
do the same with both the spatial and time derivates si- 
multaneously). For example, take the case of a flat 5D 
target space, so that there is a single physical scalar <fi. 
The unitary gauge induced metric is = rj^+d^ d v <p. 
The first higher-derivative Galileon is the cubic, coming 
from the extrinsic curvature term 



Sk 



d 4 3 



' 1 + (dej)) 2 



(17) 



Looking at the structure of the possible higher-order time 
derivatives, the only offending term is 



1 



(d# 



C C 



K 



(18) 



Expanding the denominator in powers of (dej)) 2 we 
see that every term in this expansion is of the form 
4>4> n (V '4>) 2m for some integer m and n. Integrating by 
parts, we can express each one in terms of first time 
derivatives only: 0</> n (V</>) 2m ~ ^(0 n+1 )(V0) 2m ~ 
<^ ,l+1 ^(V0) 2m . The same can be done with the higher 
Galileons and with a curved target space (see for example 
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the Hamiltonian analysis of [26, 27] in the non-relativistic 
case). 

Conclusions: There exists [15, 16, 28] a wide range of 
novel scalar field theories with interesting properties such 
as Vainshtein screening and non-renormalization theo- 
rems in common with the original Galileon models of 
[29, 30]. These properties hold out the hope that such 
models may be of use both in particle physics, and as a 
possible way to modify gravity in the infrared. However, 
coupling such fields to General Relativity in a way that 
preserves their symmetries and second order equations 
of motion seems to be impossible [19]. Instead, Galilcon- 
like scalar fields seem to most naturally couple to dRGT 
massive gravity [12]. 

The consistency of such a proposal rests on the preser- 
vation of the hard- won ghost-free structure of the dRGT 
theory. In this letter, we have shown for the first time 
that a theory of nonlinear massive gravity coupled to 
DBI scalars in such a way as to preserve the generalized 
Galileon shift symmetry and the property of having sec- 
ond order equations of motion is ghost free. Our proof is 
based on the vierbein formulation of massive gravity, in 
which the Hamiltonian analysis simplifies. Our analysis 
shows that the dRGT-DBI system provides a consistent 
framework in which models of interest to cosmology [31] 
may be developed. 
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